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Abstract 

The spectral function of ^He is calculated with the Lorentz integral transform 
method in a large energy and momentum range. The excitation spectrum of 
the residual 3N-system is fully taken into account. The obtained spectral 
function is used to calculate the quasi elastic longitudinal (e, e') response Rl 
of '^He for q = 300, 400, and 500 MeV/c. Comparison with the exact R^ 
shows a rather sizeable disagreement except in the quasi elastic peak, where 
the differences reduce to about 10% at q=500 MeV/c. It is shown as well that 
very simple momentum distribution approximations for the spectral function 
provide practically the same results for Rl as the exact spectral function. 
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Data on electromagnetic processes on nuclei can be analyzed in a very simple way with the 
help of a spectral function (SF). The approximations involved in such an analysis are few and 
transparent. There exists an extensive literature dealing with evaluations and applications 
of the SF to exclusive, semiinclusive or inclusive reactions [||. However, only for 3-body 
nuclei exact calculations of the SF have been performed 0. A complete evaluation is very 
difficult for A > 3 since it requires the knowledge of the complete set of eigenstates for 
the (A-l)-subsystem. In fact only the (A-1) ground state is often known accurately, while 
excited states, especially those belonging to the continuum, are much less under control, if 
not completely unknown. Already for ^He one finds only approximate evaluations of the SF 
0], where the final state interaction in the residual 3N-system is neglected. So the quality 
of the approximations which make use of the SF is often obscured by its poor knowledge. 

Applying the method of the Lorentz integral transform |^ one can reduce the complexity 
of the calculation of the SF considerably. In the present work we use this method to calculate 
the full SF of ^He with the semirealistic Trento (TN) potential model (central force describing 
^5*0 and ^5*1 phase shifts up to pion threshold). The result obtained is then used to evaluate 
the PWIA longitudinal (e, e') response function Rl at intermediate momenta. The resulting 
Rls are compared with the exact ones from Ref. for the same NN potential. Such a 
comparison enables us to draw conclusions about the precision of the SF ansatz in inclusive 
(e, e') scattering within a non relativistic framework. Since ^He is the lightest tightly bound 
nucleus, these results may be significant also for more complex nuclei. 

The spectral function S{k^ E) represents the joint probability of finding a particle with 
momentum k and a residual (A-l)-system with energy E. The momentum k and the energy 
E are taken with respect to the cm. and the ground state of the A-system, respectively: 

S{k,E) = —^ Y: \{^f-';ks.mj^{JoMo))\'6{E-{Ef-'-E^)). (1) 

Here and tz are the third components of the particle spin and isospin; Ef-^ and ipf-^ 
are eigenvalues and eigenstates of the (A-l)-system; Jq, Mq and Eq are the total angular 
momentum, its third component and the ground state energy of the A-system, respectively. 
There is a certain number of sum rules the SF has to fulfill: 

J dk dE S{k, E)= Jdk n{k) = 1 , (2) 



^JdEdkk'S{k,E) = {T), {E) = ^{T)-2^. (3) 

Here n{k) is the momentum distribution of the A-particle system and (T) is the mean kinetic 
energy of a particle in the ground state. The last relation in Eq. (3) is the so called Koltun 
sum rule for the mean separation energy 0. These sum rules form a set of constraints to 
test the accuracy of a calculation of S{k, E). 

In the following we will consider the proton spectral function Sp{k,E). In this case the 
first two sum rules of Eqs. (2,3) have to be modified by an additional factor Z/A on the 
right-hand sides. 

In order to express the one-body knock out cross section in terms of the SF two ap- 
proximation are required: (i) the particle interacting with the external probe is the one 
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detected in experiment, (ii) this particle does not interact with the residual (A-l)-system 
(PWIA). With these two assumptions the exclusive or semiinclusive one-body knock out 
cross sections can be written in a factorized form a ~ C AajsjSi^i — q|, £"). Here (Tat is the 
elementary cross section of the knocked out particle, kf is its momentum in the lab-system, 
q is the momentum transfer, E is the missing energy, and C is a kinematical factor. The 
so called quasi elastic (q.e.) cross section can also be written in terms of S{k,E) under the 
above assumptions. In particular the longitudinal response entering the (e, e') cross section 
reads 



RLiq,u;)c^AGl{ql) J dk, dE S,{\k, - q\, E) 6 (^u - E - ^ - ^-^^ , (4) 



where u is the energy transfer, = — cu^, Gp is the free proton electric form factor 
0, while Ma-i and k^-i = q — kf are mass and recoil momentum of the (A-l)-system, 
respectively. Here we do not consider an off-shell nucleon form factor, since our aim is a 
consistent comparison to the full Rl of Ref . , where such effects were not considered. The 
definition above includes only the proton responses of the nucleus. In principle one has also 
to consider the neutron responses, but at low and intermediate q the neutron electric form 
factor is negligible ~ 1 % at = (500 MeV/c)^). 

The SF can be calculated with the Lorentz integral transform method |Q as already 
pointed out in Ref. [§]. Let us first denote the overlap of the A-hodj bound state with the 
single nucleon plane wave 

Xp/n(k;s„Mo) = (k,s„t, = ±l/2|^o''('^oMo)) . (5) 

It represents a localized state in the subspace pertaining to the residual (A-l)-subsystem. 
Written in terms of this quantity, the proton SF 

Sp{k,E) = —^ Y: \{^Pf-'\Xr>{^■,s^,Mo))mE-{Ef-'-E^)), (6) 

looks similar to a response function of the (A-l)-subsystem with OipQ~^ replaced by Xp- 
Therefore we can proceed by analogy with the calculation of a response function. We obtain 
Sp{k, E) as a solution to the integral equation 

whose right-hand side is given by 

%ik,aR,ai) = — ^ (^p(k;s„Mo)|^'p(k;s„Mo)), (8) 

^ Sz,Mo 

where is a localized solution to the 3-body inhomogeneous equation 

{Ha-1 - E^ -an + taj)^p{k; s,, Mo) = Xp(k; s„ Mq) . (9) 

We solve Eq. expanding in hyperspherical harmonics. Complete convergence of the 
expansion is reached with similar values for the expansion parameters as in Ref. |^. The 
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inversion of the Lorentz integral transform, Eq. is carried out as described in Ref. [Q. 
Quite a good stability of the inversion results is observed. As previously we check the quality 
of the results with the help of sum rules as well. To this end we evaluate the sum rules of 
Eqs. (2,3) by an explicit integration of the properly weighted calculated SF. We obtain 
the following relative differences 0.9 % (norm), 0.2 % ((T)), and 0.8 % (Koltun sum rule 
assuming that Sn{k, E) = Sp{k, E)). Since the sum rules weight S{k, E) in different regions 
these results point out that the SF is calculated with a satisfying precision. 

Before coming to the SF, in Fig. 1 we show the n{k) of ^He for the TN potential in 
comparison to that obtained for a realistic potential (Argonne fig + Urbana IX) |^. One 
sees a rather good agreement up to almost 2 fm^^. However, different from the realistic 
result the semirealistic n{k) is considerably smaller at higher k. Most of these differences 
are presumably explained by the missing tensor force in the TN potential (see Ref. [p!0[]). 
We also show in Fig. 1 two partial momentum distributions. They are obtained from Eqs. 
(1,2) if the sum over / in Eq. (1) is restricted either to the triton ground {ntp{k)) or its 
continuum state {nt*p{k)). As expected (see Ref. |]l|) ntp{k) governs the lower k momentum 
distribution, while nt*p{k) dominates at higher k. The integration of ntp{k) leads to the so 
called spectral factor. For the TN potential one finds a spectral factor of 0.89, whereas with 
the above realistic potential a value of 0.84 is obtained. 

In Fig. 2 we show Sp{k,E). Only energies above the breakup threshold E^^ of the 
rest nucleus are illustrated, while the contribution from the bound state of the rest nucleus 
is identical to the ntp{k) of Fig. 1 (see Eq. (1)). The values of S{k,E > Ef^'^ + 1 MeV) 
are plotted in the figure. We note that S{k, E^~^) = 0, and thus S{k, E) exhibits a rather 
strong slope at low energy. For momenta below 2 fm~^ one finds a sharp maximum at about 
2 MeV above E^~^. On the contrary S{k,E) is flat in most other regions. Only for k > 2 
fm ~^ there is a ridge where the peak position shifts to higher E for increasing k. 

As already mentioned one of our aims is a comparison of the q.e. Rl with the exact 
one in an intermediate q range for the same NN potential. In Fig. 3 we show both RlS 
in comparison to experimental data. We would like to point out that the full results are a 
bit different from those in Ref. p for two reasons: (i) in the calculations of Ref. p G„(g^) 
entered erroneously with a negative sign leading to small - but not totally negligible - effects 
on Rl (e.g., peak height and high-energy tail become a bit lower); (ii) different from Ref. 
here we account for the small overbinding of the TN potential for ^He. The threshold energy 
reads uthr = -Eo(^H) — i?o(^He) + g^/2M(*^He) and we correct the overbinding by shifting 
our response to lower energies to make Uthr correspond to the one with the experimental 
Eo(^He). 

Of course these two modifications do not change the general picture given in Ref. [Q. 
The good agreement with experiment becomes even better for q.e. peak and high-energy 
tail. At low energy there is a slight improvement at g = 300 MeV/c. For the two higher q 
the agreement with experiment in the threshold region is still satisfying but not as excellent 
as shown in Ref. 

Figure 3 shows, as expected, that at g = 300 MeV/c the PWIA does not lead to a 
good description of Rl- The q.e. peak is shifted to higher energies by 15 MeV and the 
peak height is overestimated by more than 40 %. The overestimation becomes worse with 
increasing energy, while at low energy Rl is underestimated. At the two higher q the peak 
is shifted by 12 MeV, but since the peak width grows with increasing q this shift is a minor 
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effect. The shift of the peak can be quahtatively understood considering a nucleon at rest 
in a potential weU: uo can be estimated as q^/{2m) + Vf — Vi, where Vij are the potential 
energies before or after interaction with the virtual photon. While Vf is negative, it becomes 
zero in PWIA leading to an increase in u. The peak height improves with overestimations 
of 25 % at g = 400 and 13 % at 500 MeV/c. Thus one has to expect that beyond 500 MeV/c 
the PWIA is a good approximation at the q.e. peak. Beyond the peak the PWIA result still 
overestimates the exact one considerably, but the discrepancy decreases with increasing q. 
At low energy, however, the underestimation remains considerably large. 

It is advantegeous to have a simple and good approximation for the PWIA response. 
From Fig. 2 it is evident that at low k almost all the strength of S{k, E) with the disin- 
tegrated rest nucleus is found close to the breakup threshold. This suggests the following 
approximation 

Sik, E) ~ ntjXk)5{E - Eo{^K) + Eo('He)) + nt^,X^)5{E - Et^p + Eo^Re)) (10) 

for calculating inclusive processes, where Et*p is the breakup energy of the rest nucleus. One 
obtains an even simpler approximation considering that Et*p ~ Eq{^}1): 

Sik, E) ~ n{k)5{E - Eo('H) + Eo('He)) . (11) 



Eq. (11) was used e.g. in Ref. |T^ where large deviations from the full GFMC response for 
a realistic potential at g = 400 MeV/c were reported. In Fig. 4 we show the PWIA results 
with the above two approximations at g = 500 MeV/c relative to the full SF result. It is 
readily seen that the three responses are very similar, particularly in the q.e. peak region 
(at q = 300, 400, and 1000 MeV/c one has very similar results). It is worth mentioning that 
our PWIA result at g = 400 MeV/c is essentially the same as the one in Ref. |T^. This 



shows again that a semirealistic central force leads for practically to the same result as 
a realistic potential. 

In this work we obtain for the first time the full spectral function of ^He. A semirealistic 
NN potential is used. The final state interaction in the residual system is taken into account 
completely by the Lorentz integral transform method. The SF is then used to calculate the 
q.e. longitudinal response function of ^He which is compared to the exact one of Ref. 0. In 
the peak region the differences decrease with growing momentum transfer up to about 10 
% at g = 500 MeV/c, but one still finds sizeable differences apart from the peak. Similar 
results were found in Ref. [T^ for and ^He for the studied momentum transfers of 300 and 
400 MeV/c. Different from the 3-body system ^He already resembles some aspects of more 
complex nuclei and thus the general picture of the q.e. response should not change much in 
such systems. We show as well that the simple momentum distribution approximations for 
the SF provide results for Rl which are quite close to those obtained with the full SF. 

Two of us (W.L. and G.O.) thank the Institute of Nuclear Theory at the University of 
Washington for its hospitality and the Department of Energy for partial support during the 
completion of this work. The work of V.D.E. was supported by INFN and RFBR (grants 
no 96-15-96548 and 97-02-17003). 
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FIGURES 



FIG. 1. Total (solid curve) and partial momentum distributions ntp (dotted curve) and nt*p 
(dashed curve) of ^He with TN potential; also shown total result (full dots) and ntp (open squares) 
with Argonne vis + Urbana IX [9] . 

FIG. 2. Sp{k, E) of ^He with TN potential in units of fm^MeV'^ 

FIG. 3. Rl of '^He with TN potential: PWIA results according to Eq. (4) (dashed curves) and 
full results (solid curves); experimental data from Bates [11] and Saclay [12]. 

FIG. 4. Rl with SF of Eq. (10) (dashed curve) and of Eq. (11) (dotted curve) relative to Rl 
with full SF (the q.e. peak is marked by an arrow). 
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